Coherent Quantum Engineering of Free-Space Laser Cooling 
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We perform a quantitative analysis of the cooling dynamics of three-level atomic systems inter- 
acting with two distinct lasers. Employing sparse-matrix techniques, we find numerical solutions 
to the fully quantized master equation in steady state. Our method allows straightforward deter- 
mination of laser-cooling temperatures without the ambiguity often accompanied by semiclassical 
calculations, and more quickly than non-sparse techniques. Our calculations allow us to develop an 
understanding of the regimes of cooling, as well as a qualitative picture of the mechanism, related 
to the phenomenon of electromagnetically induced transparency. Effects of the induced asymmetric 
Fano-type lineshapes affect the detunings required for optimum cooling, as well as the predicted 
minimum temperatures which can be lower than the Doppler limit for either transition. 

PACS numbers: 32.80.Pj, 42.50.Vk, 32.80.Wr 



Laser cooling has led to tremendous achievements, 
enabling extremely precise measurements (see, for ex- 
ample, Ref. jDh the development of highly accurate 
atomic clocks (2J, and the production and manipula- 
tion of quantum-degenerate gases |3|. The atomic in- 
ternal structure greatly restricts the species for which 
laser cooling can be applied successfully, and determines 
whether ultracold temperatures (on the order of a few 
/iK) can be achieved. The basic Doppler cooling mech- 
anism 0, Q , based on momentum transfer from a near- 
resonant light beam, has a minimum temperature (typ- 
ically a few hundred /iK) proportional to the scattering 
rate (linewidth) of the atomic transition, and given by the 
balance between friction (cooling) and diffusion (heating) 
in the atom-laser interaction. Lower temperatures can 
be achieved either by exploiting the multilevel hypcrfinc 
structure (Sisyphus [E 0, Raman @, VSCPT Q), by 
modifying the atomic scattering (spontaneous emission) 
for example in a near-resonant cavity |lfjj |. or, more re- 
cently, by using weakly allowed, narrow two-level optical 
transitions in a second-stage for atoms that have been 
pre-cooled by other means [ll],[l2- Here, we explore cool- 
ing by dissipative radiation-pressure forces |5| in three- 
level systems under two-color excitation. Similar schemes 
have been studied previously for trapped ions 0, 0] . 
In the electromagnetically-induced-transparency (EIT) 
regime, e.g., with a weak "probe", and a strong "dress- 
ing" laser, an effective two-level system is engineered that 
allows sideband coolin g of ions to the ground state of the 
trapping potential |l4Lll5j. 

This Letter presents a general method that uses co- 
herent quantum interference to tailor an atom's internal 
energy structure that makes it more amenable to simple 
Doppler cooling in free space. To exploit the effects of 
quantum interference, at least three internal atomic lev- 
els and two lasers with distinct frequencies must be uti- 
lized. There are three basic types of three-level systems 
- the A-, V-, and S-systems — each classified accord- 



ing to the ordering of the bare quantum states in energy, 
and the possible decay pathways. The A-configuration is 
commonly used for studying of EIT and related phe- 
nomena. Here we focus on a S-type system, also known 
as a cascade configuration. We note, however, that our 
general conclusions can be applied to any type of three- 
level system, and can be extended to systems with more 
than three levels. In particular, we will focus our discus- 
sion on a cascade system suggested by Magno et al. 17'] as 
a simpler second-stage cooling scheme for alkaline-earth 
atoms, which has recently been demonstrated for mag- 
nesium |18| . The level structure and internal parame- 
ters for this system are depicted in Fig. ^ The cooling 
scheme is more effective for A-type three-level systems 
having metastable lower states, but also works for ladder 
systems with an upper state narrower than the intermedi- 
ate one ^3 ■ This cooling technique seems well suited for 
alkaline-earths, which are good candidates for the next- 
generation optical atomic clocks, studies of ultracold col- 
lisions, optical Feshbach resonances [l^, and achieve- 
ment of quantum degeneracy poj . Other elements with 
similar structure include Zn, Cd, and Hg. For lighter 
alkaline-earth atoms, such as Ca and Mg, for which tech- 
niques other than Doppler cooling are impossible or too 
difficult to be used 17], it will facilitate trapping in opti- 
cal potentials and lattices. Our treatment solves the fully 
quantum-mechanical equations of motion for the three- 
level atom in the field of both lasers. The dynamics of this 
system are determined exactly using numerical methods, 
giving quantitative predicted laser-cooling temperatures. 
The results of these calculations then provide us with 
a basis for a more intuitive interpretation of the cooling 
mechanism. Along these lines, we develop a dressed-state 
picture of the system dynamics that makes qualitative 
predictions that agree well with the numerical calcula- 
tions. 

Fig- D shows the internal atomic states in order of in- 
creasing energy as |0), |1), and |2). The transition en- 
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FIG. 1: Left side: Atomic configuration for the three-level 
cascade system in 24 Mg, utilizing the (3s 2 ) 1 5'o — * (3s3p) 
l Pi — > (3s3d) 1 i?2 transition, with levels denoted by |0), |1), 
and |2), respectively. The quantum interference between var- 
ious excitation pathways in this system can be used to in- 
crease the effectiveness of laser cooling. The energy of the 
lower and upper transition is hu^ and hu^\ respectively. 
The frequency and detuning for each laser is labeled by un 
and Si, respectively, for i = 1, 2. The spontaneous-emission 
linewidths of the states (1) and \2) are Fi and Ti, respectively. 
Right side: Characteristics of the dressed 24 Mg cascade sys- 
tem, with si(5i) = 0.001 and S2($2) = 1. Real (top graph) 
and imaginary (bottom graph) parts of the eigenvalues of the 
Hamiltonian in Eq. J3J , with dressed atomic states labeled for 
the manifold. The real parts are the energies and the imagi- 
nary parts are the effective linewidths of the dressed atomic 
system. Both are plotted as functions of the dressing laser 
detuning 82 , with fixed Si — 0. 



ergy of the lower transition, |0) — > |1), and of the upper 
transition, |1) — ► |2), are HlOq 1 ^ and Hu>q , respectively. 
We include two lasers, of frequency lo\ and oj 2 , and de- 
fine their detunings from the appropriate atomic transi- 
tions as Si = LL>i — lOq \ for i = 1,2. The intensities of 
lasers 1 and 2 are characterized by their Rabi frequencies 
fix = - (0|d|l) • Ei(x) and Q 2 = - <l]d|2> • E 2 (x), re- 
spectively, where d is the electric-dipole operator of the 
atom and Ej is the electric-field amplitude for the ith 
laser. The spontaneous-emission linewidths of states |1) 
and 1 2) are T ± and T 2 , respectively (88 MHz and 2.2 MHz, 
respectively, for 24 Mg). The time evolution for the atom 
moving in the laser field, with mass m and center-of-mass 
momentum operator p, is described by the master equa- 
tion, 

p{t)= % -\p,H] + C\p], (1) 

where p is the reduced density matrix of the atom system, 
the vacuum photon field degrees of freedom having been 
traced over, H = p 2 /2m + fiw^ |1) (1| + huffi |2) (2| + 



^iisL( x ) + ^liser ( x ) i s tue Hamiltonian of the atom- laser 
system, and the superoperator Liouvillian C describes 
effects due to coupling of the atom to the vacuum photon 
field, resulting in spontaneous emission |2lj| . 

Eq. treats all the atomic degrees of freedom (inter- 
nal and external) quantum mechanically, so its solutions 
generally provide an accurate description of the atom's 
dynamics. We thus avoid much of the ambiguity and 
difficulty associated with semiclassical approximations 
of the system. For reasonable temperatures, however, 
the number of numerical basis states required to directly 
solve the problem, even in one dimension, is impractical 
for most computers. Here we resolve this problem with- 
out resorting to Monte Carlo methods by noting that 
the matrix for the linear system equivalent to Eq. Q is 
very sparse — i.e., only a small fraction of its elements 
are nonzero. The particular structure shared by Liouvil- 
lian operators C describing relaxation processes simplify 
the numerics. Starting from the microscopic properties 
of atomic operators comprising this so-called Lindblad 
form of C, we construct a matrix with the zero elements 
eliminated. The steady-state solution of Eq. is then 
found using standard sparse-matrix inversion or propa- 
gation techniques. The result is an exact direct solution 
of a fully-quantized master equation. 

The steady-state density matrix has been calculated 
this way in one dimension. The average kinetic energy 
is then (p 2 /2m} = Tr(pp 2 /2m), which we equate with 
I^T, where T is the temperature and Ub is Boltzmann's 
constant. The large parameter space of the three-level 
atom two- laser problem has been explored, e.g., vary- 
ing the detunings, <5i and S2, as well as the strengths 
f?i and il 2 of the two lasers. However, we find that the 
laser strengths are better characterized in terms of how 
strongly they dress the atom. The relevant non-resonant 
saturation parameters, Si and s 2 , for the respective tran- 
sitions, are 

1 n 2 

Sl = H 2 + (ry2)2- (2) 

The parameter space has three distinct regimes. In the 
first, with si > 1 and s 2 arbitrary, only heating occurs, 
as expected from Doppler-cooling theory since the lower 
transition is driven strongly. In the second, with si C 1 
and s 2 <^ 1, cooling occurs only to the Doppler limit for 
the lower transition T^^ = hTi/2ks, and only in the 
range near Si — —T\/2, In this case, laser 2 has no ef- 
fect, as this amounts to simple two-level Doppler cooling 
on the lower transition with laser 1. In the final regime, 
when si < 1 and s 2 ^ 1, cooling occurs down to sub- 
stantially below Tp D ^ . Figure|21illustrates cooling in this 
regime for the particular case of 24 Mg mentioned above, 
with the temperature normalized to Tp D \ and with the 
bare two-photon resonance (Si + S 2 = 0) denoted by the 
dashed line. The steady-state temperature is plotted as 
a function of Si and S 2 for Si(<5i) = 0.001 and s 2 (S 2 ) = 1. 
Note that the saturation parameters are being held fixed 
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as the detunings are varied, so that the Rabi frequencies 
are being continuously adjusted. We see the lowest tem- 
peratures, on the order of 1CP 2 Tp D \ in the quadrant 
with 5i > and S2 < 0, as well as less extreme cooling in 
other regions. Observe that the lowest temperatures are 
obtained for frequencies detuned to the blue of the two- 
photon resonance. This seems counterintuitive, since a 
red detuning is usually required in order to have a net 
decrease of atomic kinetic energy in a photon-scattering 
event. Qualitative understanding of the cooling mecha- 
nism emerges from analysis of the simpler Hamiltonian, 
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Its complex eigenvalues have real dressed energies, and 
imaginary parts giving dressed-state linewidths — that 
is, a measure of the coupling of the dressed states to 
the photon vacuum. These dressed energies and widths 
are plotted on the right side of Fig. ^ a s functions of 
<5 2 , for the same parameters used in Fig. |21 The can- 
cellation of one of the widths can be viewed as an EIT 
effect: since laser 1 is perturbative while laser 2 strongly 
dresses the upper transition, the new eigenstates of the 
system, denoted |+) and |— ), are well approximated as 
linear combinations of the bare states |1) and |2). These 
states have the modified energies and widths shown in 
Fig. H The linewidth modifications can be viewed as a 
Fano interference |23| , in which the dressing laser transi- 
tions caused by the probe laser enable multiple coherent 
pathways among the bare states. Constructive or de- 
structive interference respectively increases or decreases 
the atomic linewidth. 

The cooling mechanism is thus qualitatively explained 
as ordinary two-level Doppler cooling. But instead of us- 
ing a transition between two bare states of an atom, the 
transition occurs between a (mostly) unmodified ground 
state, and a dressed excited state, with a shifted energy 
and a new linewidth that can be narrower than the bare 
linewidth of the lower transition. As the probe laser is 
scanned, the detuning relative to the dressed energy lev- 
els is varied, but since these levels are shifted from their 
bare energies, resonance occurs for different detunings 
than are encountered in the bare system. In fact, the 
shifts of the eigen-energies in Fig. ^ from the bare ener- 
gies explain the apparent observation of blue two-photon 
cooling in Fig. [21 In the dressed system, the bare two- 
photon resonance is no longer meaningful, and the cool- 
ing region is in fact to the red of a dressed resonance. 

An additional caveat applies when mapping this sys- 
tem onto Doppler cooling theory: the lineshapes are not 
Lorentzian, but are in fact asymmetric Fano lineshapes 
for the dressed system, as shown in the upper part of 
Fig.[3]as a function of 5i for a fixed value of 82 = —T\/2. 
This fact changes the optimum-detuning condition into a 
new one: maximum cooling for a given transition occurs 
when the probe laser is detuned from the dressed excited 
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FIG. 2: (Color online) Steady-state laser-cooling tempera- 
tures for a 24 Mg three- level cascade system, as a function of 
the two detunings <5i and S2, obtained from exact numeri- 
cal solutions of Eq. Q. The detuning-dependent saturation 
parameter for the lower transition (perturbative probe laser) 
is si(<5i) = 0.001 for both plots, and for the upper transi- 
tion (dressing pump laser) is 52(^2) = 1 and 5 in the upper 
and lower plot, respectively . The temperature is normalized 
to the one- dimensional Doppler limit for the lower transition 
T^ 113 ' = 7hTi/A0kB p^l. which is the optimum temperature 
expected for cooling with just one laser. The dashed line in- 
dicates the location of the bare two-photon resonance. This 
example illustrates the main parameter regime where sub- 
Doppler cooling occurs. Note that, in order to hold the sat- 
uration parameters fixed as the detuning is varied, the Rabi 
frequencies for i — 1,2, are continuously adjusted. 
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FIG. 3: Upper plot: Absorption spectrum (solid line) as a 
function of 5i , illustrating the asymmetric lineshapes for the 
dressed system, for a fixed dressing-laser detuning b~2 = —Ti. 
The peak of each lineshape is located at a dressed eigen- 
energy, and lineshapes at the same energies and widths, but 
with Lorentzian (symmetric) linshapes are plotted with dot- 
ted lines. As an example, optimum-cooling detunings relative 
to the leftmost resonance are illustrated with a black arrow 
for the true asymmetric lineshape and with a gray arrow for 
the hypothetical symmetric-lineshape case. The value of the 
optimum detuning, as well as the slope of the lineshape is seen 
to be different for these two cases. Lower plot, solid line: the 
ratio of the maximum slope of a Lorentzian lineshape with 
width Fi to the slope of the asymmetric lineshape, as a func- 
tion of 61 with 82 = —T\/2. This ratio provides an indication 
of the expected cooling for the dressed system relative to the 
Doppler limit for the lower transition. For comparison, fully 
quantum numerical results are indicated by data points. 



state precisely to the inflection point of the absorption 
spectrum. This can be understood by noting that the 
force / applied to the atom due to the laser beam is pro- 
portional to the absorption rate, for a given 5±. As is 
often utilized in semiclassical cooling theories, the fric- 
tion coefficient a for the atom moving in the laser field 
is given by 



dv 



(4) 



where v is the atomic velocity. Since the detuning of 
the laser and the resonant atomic velocity are linearly 
related, the derivative of the absorption spectrum with 
respect to S± also yields a maximum in the cooling force. 
This is evident in normal Doppler cooling because the 
optimum detuning occurs when S = — F/2, which is the 
inflection point of the Lorentzian. In general then, for 
asymmetric lineshapes, the optimum detuning does not 
obey such a simple relation, but depends on the degree 
of asymmetry. 

From this complete picture of the cooling mechanism 
as a weak probe applied to a dressed three-level system, 
the minimum expected temperatures can now be deter- 
mined, allowing for the detuning modification due to 
asymmetric lineshapes. The lower part of Fig. 02 the ra- 
tio of the maximum slope of a Lorentzian lineshape with 
width Ti to the slope of the asymmetric lineshape, as a 
function of Si with 62 = —T\/2. For comparison, fully 
quantum numerical results are indicated by data points. 
This ratio provides an indication of the expected cooling 
for the dressed system relative to the Doppler limit for 
the lower transition. Note that the expected tempera- 
ture, due to the asymmetric lineshape, is predicted to be 
lower than the upper-transition Doppler limit, indicated 
by the horizontal dotted line in the lower plot of Fig. [3] 
This prediction is supported by the numerical data. 

In conclusion, coherent engineering of an atomic three- 
level system can optimize the effectiveness of two-level 
Doppler cooling. By strongly driving a particular tran- 
sition between two excited internal states of the atom, 
dressed energy eigenstates are created with modified 
linewidths which, due to mixing, can vary anywhere in 
the range between the smallest and the largest of the two 
bare linewidths. Smaller linewidths lead to lower temper- 
atures, but the additional effect of asymmetries in the 
linewidths of the dressed states can lead to even lower 
temperatures, below the Doppler limit of the upper tran- 
sition. The ability to tailor the degree of cooling lends 
this technique additional utility. A dressing scheme can 
be suited to the characteristics of a particular atom, and 
real-time adjustment of the cooling properties can allow 
narrowing of the velocity-capture range as an atomic gas 
is cooled. Utilizing such coherent effects should lead to 
relatively simple schemes for cooling far below the typical 
Doppler limit. 
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